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We give a detailed study of the asymptotic behavior of field commutators for linearly polarized, 
cylindrically symmetric gravitational waves in different physically relevant regimes. We also discuss 
the necessary mathematical tools to carry out our analysis. Field commutators are used here to 
analyze microcausality, in particular the smearing of light cones owing to quantum effects. We 
discuss in detail several issues related to the semiclassical limit of quantum gravity, in the simplified 
setting of the cylindrical symmetry reduction considered here. We show, for example, that the small 
, G behavior is not uniform in the sense that its functional form depends on the causal relationship 

between spacetime points. We consider several physical issues relevant for this type of models such 
as the emergence of large gravitational effects. 
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I. INTRODUCTION 

> 

in 

' Among the many symmetry reductions of General Relativity that have been considered in the past, linearly polarized 
cylindrical waves (also known as Einstein- Rosen waves 0,0) have been the focus of intensive study. They provide a 
model with an infinite number of degrees of freedom that can be exactly solved in spite of the fact that it is non-linear. 
This is not only true classically but also quantum mechanically, and hence this system is a valuable tool to explore 
the physics that may be found if a successful quantization of full general relativity is ever achieved 0, 0, El • 

One of the main reasons behind this success is the fact that the physical Hamiltonian is a function of the free 
Hamiltonian of a 2+1 dimensional, axially symmetric, massless scalar field evolving in an auxiliary Minkowskian 
background ID, • In a previous paper 0] we took advantage of this fact to study the quantum corrections to 
5h i the spacetime structure by considering the commutator of this scalar field at different spacetime points. As is well 
known, microcausality in quantum field theories can be discussed by looking at the commutator of quantum fields (or 
anticommutator in the case of fermions) . In the standard examples the microcausality requirement means that this 
(anti) commutator must vanish for spatially separated spacetime points. A similar argument can be made for vector 
fields, though issues of gauge invariance change some of the conclusions; in particular if one computes the commutator 
of the four-vector potential at two spatially separated points, it may be different from zero in some gauges, even 
though it is always true that the commutator of gauge invariant objects is zero for such points. 

The issue of gauge invariance in the context of cylindrical gravitational waves has been discussed at length by Bicak 
and collaborators ■ These authors show that is is legitimate to use the Ashtekar-Pierri gauge fixed action written 
in terms of the scalar field that encodes the physical information in this model to obtain gauge invariant structures 
such as Dirac observables or the S matrix. This justifies the computation of the type of objects -field commutators- 
that we will be considering here to extract conclusions about the quantum structure of spacetime. 

The particular problem that we will be concerned with in this paper is the detailed study of the field commutator 
(given as a certain integral) and, in particular, its limiting behavior when the time and length separations are much 
larger than the natural length scale of the problem the Planck length-. In order to do this the procedure of expanding 
the integrand of the field commutator as a power series in the gravitational constant and other asymptotic parameters 
is not useful. In fact it will be necessary to adapt some methods developed for the asymptotic analysis of integrals 
and get a consistent procedure to expand the relevant objects as asymptotic series. It is possible to understand this 
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a posteriori as a consequence of the fact that some limiting behaviors (i.e. in G) of the field commutator change in 
a non trivial way as the spacetime intervals go from spacclikc to timclikc or when one of the spacetimc points lies in 
the symmetry axis. Also the functional dependence in some of these parameters is highly non-polynomial. This is 
not what one would expect to obtain in the familiar perturbative treatment of quantum field theories (QFT's). 

The paper is divided in two main sections: a physical discussion of the behavior of the field commutator followed 
by a detailed description of the asymptotic methods necessary to study the different physically relevant regimes. 
Specifically, after this introduction we will give the different asymptotic expansions for all the relevant parameters 
(involving the G — > limit and also the limits in which the difference in the time coordinates or the radial coordinates 
go to infinity) . Using them we will discuss the main physical consequences of the quantization of this model as far as 
microcausality is concerned. In this respect it is particulary interesting to point out the existence of a certain type of 
large quantum mechanical effects (in a sense that will be made precise later) when one of the spacetime points in the 
field commutator lies in the symmetry axis. 

A technical issue that should be considered is the role of regulators in the final physical results. As is well known 
regulators are necessary to give sense to otherwise ill-defined objects. They must be introduced, for example, to obtain 
a finite norm vector by acting with the field operator on the vacuum state. Cut-offs are a simple way to regulate 
amplitudes. It is conceivable that physical regulators exist that restrict, for example, the integration intervals for some 
physical objects (such as the field commutators considered here) to finite real intervals. However it is also possible 
that they are just a convenient way to render some physical objects finite in such a way that no footprint is left in 
the final results. This is the philosophy that one has in mind in the usual renormalization scheme where cut-offs are 
taken to infinity and disappear from physical quantities. 

The point of view of this paper is that the un-regulated objects (integrals), when defined, are a good approximation 
to the regulated ones. The asymptotic analysis of regulated commutators and their relation with the un-regulated 
ones will be discussed elsewhere. 

In the usual perturbative QFT computations Green functions, matrix elements, and similar objects are expanded 
as power series in the coupling constants of the model with coefficients that are usually written as regulated integrals. 
This is a necessary step because it is usually not possible to write them in closed form. Here the situation is different 
because it is possible to write down expressions for the objects of interest (field commutators in the present example) 
that depend on the coupling constants in a non trivial way. This has the advantage of allowing us to use approximation 
techniques specially adapted to their specific form and much better suited to its study. It also permits to consider 
some problems that may be difficult to tackle for the usual QFT's; for example one can try to figure out if the 
asymptotic behavior of some regulated object coincides with its asymptotic behavior after the regulator is removed. 

A detailed technical discussion of the techniques needed to explore the relevant asymptotic regimes is presented in 
section ITTT1 These techniques will be useful tools in order to analyze other types of physical objects (such as S-matrix 
elements) so they provide the necessary background for a perturbative framework properly adapted to this model, this 
is why we invest some time to study them here. We end the paper with a summary of the main results, comments, 
perspectives for future work, and several appendices. Numerical and algebraic computations have been performed 
with the aid of Mathematical . 



II. THE FIELD COMMUTATOR 



We start with a brief description of the model and introduce our conventions and notation. As is well known 
QUI Einstein- Rosen waves correspond to topologically trivial spacetimes with two linearly independent, commuting, 
spacelike, and hypersurface orthogonal Killing vector fields. The metric in this case can be writen as 

ds 2 = e^-^{-alT 2 + alR 2 ) + e -^R 2 d9 2 + e^dZ 2 , (1) 

where we are using coordinates (T, R,9, Z), T E R, R E [0, oo), 6 E [0, 2ir), Z € K, and ip and 7 are functions only 
of R and T. The Einstein field equations for this metric are very simple: the scalar field ip must satisfy the wave 
equation for a massless, axially symmetric scalar field in three dimensions 

d 2 i> - d 2 ^ - ^8bM> = 

and the function 7 can be expressed in terms of this field on the classical solutions 0, E| as 

l(R) = \j*dRR [(d T V) 2 + (d a if>) 2 ] . 
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Throughout the paper we will use a system of units such that c = K = 1 and define G = /1G3, where G3 denotes 
the gravitational constant per unit length in the direction of the symmetry axis 1 . The metric function j(R) admits 
a simple physical interpretation: apart from a factor of 8G it is the energy of the scalar field in a ball of radius R, 
and 7oo = liniR^oo j(R) the energy of the whole two-dimensional flat space. Furthermore, 700 /8G coincides with the 
Hamiltonian Hq of the system obtained by linearizing the metric Q 3, |6j . 

In order to arrive at a unit asymptotic timelikc Killing vector that allows us to introduce a physical notion of energy 
(per unit length) it is convenient to use coordinates (t,R,9,Z) defined by T = e >°°/ 2 t. In these coordinates the 
metric takes the form |2|, |9( 

ds 2 = e^{-e^dt 2 + dR 2 ) + e^R 2 d9 2 + e^dZ 2 . 

By choosing a metric function ijj with a sufficiently fast fall-off as R — > 00, this metric describes asymptotically flat 
spacetimes with a certain (non-zero) deficit angle and such that d t is a unit timelike Killing vector in the asymptotic 
region. 

The Einstein field equations can be obtained from a Hamiltonian action principle 0, [ToL 1 1 1 1] where the Hamiltonian 
H is a function of that corresponding to the free scalar field, H n : 

ff = £(if„) = i(l- e - 4Gfl »). 

In the following we will refer to t as the physical time and to H as the physical Hamiltonian. 

In terms of the T-time and imposing that ip be regular at R = the classical solutions for the field i/j can be 
written as 

/>oo 

iP(R,T) = V4G / dkJ (Rk)[A(k)e- ikT + A^(k)e ikT ] 
Jo 

where A{k) [and its complex conjugate AT(fc)] are determined by the initial conditions 2 . The free Hamiltonian Hq 
can be written now as 

/>oo 

7oo = H = / dkkA^i^Aik). 
Jo 

Using this expression, we can obtain the t-evolution of the field 

POO 

4> E (R,t) = V4G dkJo(Rk) A{k)e- lkte ~''°° +A\k)t 



ikte 



The quantization can be carried out by following the usual steps. We introduce a Fock space in which ^(R, 0), the 
quantum counterpart of ijj(R,0), is an operator- valued distribution Its action is determined by those of A(k) 

and A^ (k) , the annihilation and creation operators with non- vanishing commutators given by 



A(k x ),AHk 2 ) =s(k 1 ,k 2 ). 



Explicitly, 



Tp(R,0)=4> E (.R,0) = V4n / dkJo(Rk) i(fc) + i f (fc) 

Jo L 



Evolution in T is given by the unitary operator Uq(T) = exp(— iTHo) where 



Ho 



dkkA\k)A(k) 



is the quantum Hamiltonian operator of a three dimensional, axially symmetric scalar field. The quantum scalar field 
in the Heisenberg picture is hence given by 



ij{R,T) = Ul(T)i>{R,0)Uo{T) 



4G I dkJo(Rk) A{k)e- lkT + A\k)t 
'0 



1 Notice that in this system of units G has dimensions of length. 

2 Throughout the paper improper integrals are understood in the Riemann sense = lim^— tao J ■ 
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If we choose the physical time t to define the evolution in our model, the quantum Hamiltonian is H — E(Hq) = 
— e~ °) and unitary evolution is given by U(t) — exp(-itH). With this time evolution the annihilation and 
creation operators in the Heisenberg picture are 



A E (k,t) = tj\t)A(k)U{t) 



exp 



-itE(k)e 



-iGH 



A(k), 



A%(k,t) = A\k) exp itE(k)e 



-AGH 



where E(k) — — e 4Gfc ), and the field operator evolved with the physical Hamiltonian [that we denote as $ E (R, t)] 
is given by 

4> E {R,t) =ViG / dkJo(Rk) A E (k,t) + A E {k,t) 



The field commutator \4>e(Ri, ti), ^e{R2i ^2)] can be computed from these expressions One of its interesting 
features is the fact that it is not proportional to the identity as in the case of free theories (or if we consider the 
quantum evolution defined by the Hamiltonian Hq of the linearized model) but is a non-diagonal operator in the 
chosen basis. We are hence led to consider its matrix elements. We will concentrate here on the most relevant of 
these elements (at least from the perspective of the microcausality of the classical background of the model), namely 
the vacuum expectation value, which is explicitly given by 



1 

8iG { 



bpE(Ri,ti),Tp E {R2,t2)]\0) = dkJ {R 1 k)J {R 2 k) sin 



t 2 -u 

AG 



(2) 



As we can see it depends on the time coordinates through their difference t 2 — t\ (which we will assume to be, e.g., 
positive) and depends symmetrically on R\ and R 2 . The functional dependence in G is less trivial, a fact that will 
require especial attention when studying the limit in which the relevant lengths and time differences are much larger 
than the Planck length. For comparison purposes we remember that 



^(0\^(R llTl )J(R 2 ,T 2 )}\0) 



dk Jo{R 1 k)J (R 2 k) sin [k{T 2 - Ti 



(3) 



In the following it will be convenient to refer the dimensional parameters of the integral to another length scale 
that we choose as R\. We hence define R 2 = pRi, t 2 —t\ = tR\, A = Ri/AG and rewrite as 



1 A f f°° 

— (0| [i> E (Ruti),^E(R 2 ,t 2 )} 10) - ^\ J a dqJ {\q)J { P \qy TX{l - e q) 

after introducing the new variable k = q/AG. Here 9 denotes the imaginary part. 



(4) 



A. Asymptotic behavior in p: smearing of light cones 



The integral (J1J can be written as a standard ^.-transform with asymptotic parameter p, by the change of 
variables t = q\, 



J£*lJ dtJ (pt)J (t)e 



rrA(l-e-'/ A ) 



(5) 



In this case a straightforward Mellin transform analysis [1J| (discussed in section IHIf) gives the following asymptotic 
behavior in the p — ► 00 limit 
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(6) 



The asymptotic behavior in the p — > 00 limit gives a precise and quantitative description of the smearing of (cylindri- 
cally symmetric) light cones because it shows, for example, that for fixed values of Ri and (t 2 — ti) the expectation 
value of the commutator (over 8iG) falls off for large values of R 2 as 



2G{t 2 -h) 



Rl 



O 



R 5 , 
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This means that even for large spatial separations the scalar field that encodes the gravitational field does not commute 
with itself. 

In the p —* + limit, on the other hand, we get 



Ri 



dtJ Q {t)e lTX{1 - 



O(p). 



(7) 



This shows that the commutator is a continuous function of p in p = 0, a fact that will be important in the analysis 
of the semiclassical limit. 



B. Asymptotic behavior in r: large quantum gravity effects 

The integral in JSJ has the convenient form of a /i-transform if the asymptotic parameter is chosen to be p; however 
this is no longer true if the asymptotic parameter is r (which corresponds to considering large separations in the time 
coordinates). This introduces some mathematical difficulties in the asymptotic analysis that will be discussed later. 

In this case one has to perform separate analyses for p — and p ^ 0. In p = one finds that the asymptotic 
behavior when t — * oo is given by 



" A " 


1/2 , 


i 


2tt log r 





1 

R~i 

whereas for p ^ we find 
1 



f|e i ^ +rA ~ Alog '' rA ^e^ A r(iA) + e _i ^ ~ TA ~ Alo s( rA )] e - 



2-THA) } +0 



(logr) 3 / 2 



(8) 



2nRiy/p\ogT 

.y<( e i[ % + tA ~ a (!+p) lo s( rA )] e ^ A(1+p) r[iA(I + p)} + e~ i[ % - tA - a ( 1 +/°) i°sM)1 



-A(l-p)log(rA)] e f A (!-p)r[iA(I — p)] + e i [' rA_A (P _1 ) lo s('< 



A )h 



x( - 1+p) T[-i\(l + p)} 
x ^T[iX{p-l)]\ 



(9) 



O 



(logr) 2 



The most striking feature of these expressions is the unusual dependence on the asymptotic parameter r; in fact the 
dependence on inverse powers of logarithms (especially on the inverse square root of logr) cannot be obtained by 
direct application of the usual asymptotic expressions derived by Mellin transform techniques. It is also remarkable 
how slowly the commutator decays in r -in particular in the axis p = 0- a fact that is suggestive of the large quantum 
gravity effects discussed by Ashtekar in Outside the axis the decay is quicker but still rather slow. A consequence 
of the different asymptotic behaviors in r for p — and p ^ is the impossibility to recover JSJ) as the limit when 
p — * of (0 • It is interesting to point out that the frequency of the oscillations of the commutator in the r parameter 
is controlled by the value of A (proportional to the inverse of G) in such a way that although the amplitude of the 
oscillations decays very slowly they will average to the value of the free commutator on scales much larger than the 
Planck length. 

The t — > limit is simpler to analyze. Actually we find that the series obtained by expanding e lrA ( 1 - e t/A ) as 
a power series in e -t / A , exchanging integration and infinite sum, and computing the resulting integrals gives an 
expansion that converges to the value of the commutator. 



C. Asymptotic behavior in A: the semiclassical limit 



The possibility of studying the asymptotic behaviors in p and r by using the powerful mathematical tools provided 
by the theory of Mellin transforms relies on the fact that these integrals can be written as /i-transforms of the type 
mentioned above. However this is no longer possible if one wants to study the limit A — > oo (corresponding to the 
situation when R\ is much larger than the Planck length) because of the particular dependence of these integrals on 
A. There is, however, a way out of this if one is willing to abandon Mellin transforms: writing the integral in Q as a 
multiple integral by introducing the integral representation of the Bessel functions J n (n = 0, 1, . . .), 

Jn{z) = I JL e ic-D, 

V ' 2ni / 7 t n + l 
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FIG. 1: Regions in the (p, r) plane used in the discussion of the A asymptotics. Region I corresponds to < r < \p— 1|, region 
II to \p — 1| < r < p + 1, and region III to p + 1 < r. 



where 7 is a closed, positively oriented, simple path in the complex plane surrounding the origin. By doing this the 
right hand side (r.h.s.) of @ can be rewritten in the following form 



wM^fr fW dtj * 9 J-e4*<*-M<«--*>-<«-« ]}. 

#1 I 4tt 2 7 jL, J 72 tit 2 I 



(10) 



Although the contours 71 and 72 may be different in principle we will take them equal in practice. Besides, we will 
see in the next section that it is convenient to choose them in a specific way. 

As in the r case one has to perform separate asymptotic analyses for p = and p 7^ 0. In order to proceed it is 
necessary to divide the (p, r) plane into the same regions that appear in the discussion of the "free" commutator © 
which corresponds to the evolution dictated by the Hamiltonian Hq [6j (see fig. ^1 . 

In p = one finds that the asymptotic behavior when A — * 00 is given by 
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i?! 2A(l-r 2 ) 5 / 2 

If p ^ the asymptotic behavior in the different regions shown in fig. ^is the following: 
Regions IA and IB 



t J 2[1 + p 4 + 2t 2 - 3t 4 + 2p 2 (r 2 - 1)^(1 + pf - r 2 
2ttEiA \ (1 + p - t) 2 (1 - p + r) 2 (-l + p + t) 2 (1 + p + r) 2 

2r 2 



4p 



[ p4 + (t 2 _ 1)2 _ 2(1 + r 2 )p 2 ]V(l + p) 2 



4p 



(1 + P) 2 - 



(1 + P) 2 - 



1 

A 2 



Region II 



nR ly /p 



K 



(P-1) 2 



4p 



1 ( 



-if e^Ir+lp-iKi+iog^-)] 
V2vrAp|l-p|log^ 



A 3 / 2 



(11) 
(12) 



(13) 



(14) 
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Region III 




(i-p) 2 ; 

| e iA[r-|p-l|(l+lo gT ^)] e iA[ T+ ( p+1 )( log l±£_l)] 
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}+0 


_A 3 /2_ 



V27rAHl-p|log^| V27rAp(l + p)lo gT ^ 

Here K(k) and -B(fc) denote the complete elliptic integrals of the first and second kind respectively, defined by 

d9 



(15) 



K(k)= / - E(k)= / d6Vl-k 2 sin 2 9. 

Jo VI - k 2 sin 2 6> 

Some comments are in order now. The first is that the A independent terms in the above expressions correspond 
to the commutator obtained from the free Hamiltonian Hq both in the axis and outside the axis. The remaining 
terms (except when p = and r > 1) are corrections to this free commutator that fall off to zero as A — > oo, and 
have an additional, non-polynomial dependence in 1/A. Since the free commutator defines a characteristic light cone 
structure these terms are responsible for the smearing of the light cones in this model. It is worthwhile pointing out 
that the asymptotic behavior in A is different in regions I, II, III, and in the axis -this is the reason that explains 
the appearance of singularities in the borders between adjacent regions 3 - and it has a non-polynomial dependence 
on 1/A. This kind of behavior cannot possibly appear in an ordinary perturbative QFT where the relevant objects 
(propagators, Green's functions, and so on) are expanded as power series in the coupling constants. One of the novel 
features of the approach that we follow in this paper is that by using more sophisticated approximation techniques, 
and taking advantage of the fact that we have closed explicit expressions for the objects of interest, we are able to 
extract such non trivial behaviors. In the axis we notice that for t > 1 there is a A-independent contribution that 
corresponds to the free commutator and an oscillating contribution with a frequency that depends on A (this term 
is similar to the one obtained in the asymptotic analysis for r). For large values of A this term oscillates very fast 
and averages to zero. For r < 1 we get a correction to the free commutator (which is zero in region I) that goes to 
zero as 1/A. Outside the axis we see that the corrections to the free commutator fall-off as 1/A in region I -outside 
the light cone-, but only as l/VA (with the additional non-polynomial terms) in regions II and III. The presence of 
two oscillating terms in region III produces some interference effects that manifest themselves as a checkered pattern 
in plots of the commutator -especially close to the axis- which suggests a division of spacetime into cells of a size 
governed by the value of the gravitational constant (through A). These are shown in fig. |2] where we plot the free 
commutator over 8iG plus the first asymptotic correction, given by the above expressions in each of the regions. For 
comparison we also display a plot constructed with the power series representation obtained in |(| : 

v ' ^ n— v 7 71—0 v ' ) 

, . n 2 + A 2 (l+p 2 ) 
Gn{p) = 2p~X 2 

and Q_i(x) = ttF (|, |; 1; ^) /\/2x [for x > 1] is the associated Legendre function of the second kind [l^. As can 
be seen, the result of the asymptotic approximation matches that obtained with the power series expansion (|16f) 4 . In 
figs. |21 and 0] we also plot the commutator as a function of t for p — and a value of p different from zero. It is 
interesting to point out that, as r grows, a distinctive beating pattern appears due to the interference of terms in the 
asymptotic expressions for region III mentioned above. 

Comparison with the results of numerical integration confirm the accuracy of the approximation provided by the 
asymptotic expressions, as long as one is far enough from the boundary between regions. This is shown in figs. |3| 
and 01 where we compare the asymptotic approximation with the result obtained by numerically computing the field 
commutator at p = and p = 1.5. 

The limit A — > is directly obtained by truncating l|16f) to the desired number of terms. 



where 



3 As shown in |(| the only real singularity that appears now corresponds to p = 1. 

4 The very slow convergence of this oscillating series for large values of r or large values of A makes it impractical for numerical compu- 
tations. This is the reason why we only give the lower portion of the plot in fig. |5] 
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G = 0.01 G = 0.02 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 



P P 

FIG. 2: Density plot of the asymptotic approximation in A for the field commutator (over 8iG) as a function of (p, r) for 
G = 0.01 and G = 0.02. Notice the singular behavior of the approximation in the boundaries between the regions I, II, III, and 
on the axis. For comparison we show part of the plot for G = 0.02 obtained with the series expansion 116H . Sections of these 
plots for p = and p — 1.5 with G = 0.02 are shown in figs. [3]and|l] 

III. ASYMPTOTIC ANALYSIS: MATHEMATICAL DETAILS 

We discuss in detail here the techniques used to obtain the asymptotic expansions discussed in the previous section. 
Even though we will be mostly employing standard techniques it is nonetheless necessary to adapt them to the different 
parameters that appear in the problem, (p, r, A). Before studying the asymptotic expansions we present some results 
on Mellin transforms and explain the basics of their use to obtain the asymptotic behavior of integrals. 

Whenever possible it will prove useful to write the integrals under consideration as /i-transforms, i.e. objects of the 
type 

/'OC 

H[f;C] ■= / h((t)f(t)dt 
Jo 

where £ is the asymptotic parameter (that in practice approaches either or oo) and / and h are locally integrable 
functions in (0, oo) [functions that are integrable in every closed interval in (0, oo)]. This will allow us, in many cases, 
to employ standard asymptotic techniques (see e.g. |l4j ) based on Mellin transforms 5 that can be summarized in the 
following 

theorem: Let h : (0, oo) t— > C and f : (0, oo) i— ► C be locally integrable functions. Let us assume that they respectively 
have the following asymptotic expansions when t — > oo and t — > + : 

oo N(m) 

h(t) ~ e- dt " E Cmnt- rm (l0gt) n (17) 

m— n— 



5 The Mellin transform of a locally integrable function / in (0, oo) is given by J °° dt t z 1 f(t). It is a holomorphic function defined on 
the strip a < 5R(z) < (3 of the complex plane where the integral is absolutely convergent. 
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FIG. 3: Asymptotic approximation in A for the field commutator (over 8iG) as a function of r for p — and G = 0.02. We 
plot, for comparison, the commutator corresponding to the free Hamiltonian discussed in @. The solid line corresponds to a 
numerical computation of Q. As can be seen the asymptotic approximation is very good for most of the values of r. 



> 0, v > 0, 5i(r m ) grows monotonically to oo, and N(m) > is a finite integer for each value of to] and 



N(m) 



/(t)~e-«^ £ E Pmnt a -(\0gtr 



(18) 



m— n— 



> 0, /i > 0, 3?(a m ) grows monotonically to oo and N(m) > is a finite integer for each value of m]. 
Let M[h; z] and M[f;z] be the Mellin transforms of h and f (in the generalized sense, see XlA l for details), holo- 
morphic in the strips a < 5R(z) < (3 and 7 < < 5 that we assume overlapping. Let us suppose also that the 

Parseval identity holds: 6 



1 />r-Hoo 1 pr-\-zoo 

H[f; C] = 7T- / dzC z M[h; z]M[f; l- z } = — dz( 

Z7TI Jr—ioo Z7T1 J r —ioo 



r+zoo 



z G(z) 



(19) 



with r in the intersection of the holomorphicity strips. Lf Ixaiiyi^^ G(x + iy) = for each x £ [r, R], (R € WL) and 
XT' dy\G(R + iy)\ is finite then 

H[f; C] ~ - ™*(C z M[h; z]M[f; 1 - z}) 

r<9t(z)<R 

is a finite asymptotic expansion of H[f;(] as ^ — * 00 with respect to the asymptotic sequence $j,m(C) = 
{C ~ Qj (log C)" J ~ m }; m — 0j 1j ■ ■ ■ i n j 3 = 0, 1, . . . with error 0((~ R ). If the previous hypotheses hold for arbitrary 
values of R we get an asymptotic expansion with an infinite number of terms. 



Conditions for this identity to hold can be found in Il4l . 
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FIG. 4: Asymptotic approximation in A for the field commutator (over 8iG) as a function of r for p = 1.5 and G = 0.02. We 
plot, as before, the contribution corresponding to the free Hamiltonian and a numerical computation of the commutator. 



In the following we will only need to consider the case when d ^ and q = 0. The result of the previous theorem 
can be written as 



JV(m) 



m=0 n=0 j=0 

with M^[h; z] denoting the k th derivative of the Mellin transform with respect to to z. 



z—l+a^ 



(20) 



A. Asymptotic behavior in p 



We write the integral that gives the field commutator in the form JSJ and choose p as the asymptotic parameter. 
In order to apply the previous theorem we need the asymptotic expansion of f(t) = Jo(i)e lrA ^ 1_e * /A ) as t — * + and 
the Mellin transform of h(t) = Jo(t). From the first we get at = k, (k = 0, 1, . . .), poo = 1) Pio — i T , etc. The Mellin 
transform of h is given by 



M[h;z] = 



2 z " 1 r(z/2) 

r(i-^/2) 



Since M[h; k] = for even values of fc, we get the asymptotic expansion in odd powers of 1/ p given in 10. 

In order to obtain the asymptotic expansion of the commutator when p — > + we change variables in (J5J according 
to s = pt, 
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take 1/p as the asymptotic parameter, and switch the roles of / and h 0] [/(s) = Jo(s), h(s) = Jo(s)e lrA ' 1 e 1 )]. 
The Mellin transform of h is given now by 

/>oo 

M[h;z} = / dtt*- 1 J (i)e* rA(1-e 5 , 
Jo 

which converges in the strip < 3t(z) < 3/2. We do not know any closed expression for this integral but it suffices 
to obtain the first term of the asymptotic expansion of the commutator when p — > + . The asymptotic expansion 
of / gives now = 2k, (k = 0, 1, . . .), poo = 1, Pio — — 1/4, etc. and using the theorem we obtain the expansion 
(Q). Under mild restrictions 0] it can be shown that M[h; z] admits an analytic continuation for $t(z) — > oo which is 
meromorphic, at worst, in the complex plane, and that can be used to get further terms of the analytic expansion as 
p — > + . Notice, however, that beyond the first term we cannot use the integral representation of the Mellin transform 
of h given above because the integral becomes divergent for z — 2k; k > 1. This prevents us from getting further 
terms in the asymptotic expansion as p — * + by using this method. 



B. Asymptotic behavior in r 

The integral in (0J is not an /i-transform when r is chosen as the asymptotic parameter but can be written as such 
by using the change of variables u = Ae~'/ A . This gives 

\ C P ^ — IT'LL *\ 

Even though this integral has now the appropriate form with f(u) — ^X[o,A]( u )>^o(Alog v) Jo(p^\og j) and h(u) = 
e~ lu , the asymptotic behavior of /(it) as it — > + is not of the type considered in the theorem. This has some 
unpleasant consequences, the Mellin transform of / cannot be analytically continued as a meromorphic function on 
C, and we cannot use the expressions given in the theorem for the asymptotic behavior of /i-transforms. This forces 
us to follow a more complicated approach. We will consider the cases p = and p ^ separately. 



p = 



We have to study the integral 



A „-ir« , u 



when t — > oo. Let us set h(u) = e lTU and f(u) = X[o.x]( u )^Jq{^^°S j) [remember that /i-transforms are defined as 
integrals over (0,oo)]. The Mellin transforms of these functions are 



M[h;z] = T- z r(z)e-™ z/2 , 

convergent in the strip < 3?(z) < 1, and 



\z-l 

M[f;z] = 



v/A 2 + (z-l) 2 ' 

which converges in 1 < $t(z). Several comments are in order. First we see that the Mellin transform of / cannot be 
analytically continued as a meromorphic function over the complex plane because of the cuts coming from the square 
root. If we choose the branch 7 given by log(z) = log |^| + iarg(z) with arg(z) S (— tt,tt] the cuts are those parts of 
the imaginary axis with \z\ > A. Another comment is that the regions where the Mellin transforms converge do not 
overlap; this fact precludes us from directly employing the Mellin-Parseval formula (|19fl . The reason can be traced 



7 We will use this choice throughout the paper. 
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back to the behavior of h(t) as t — ► + ; fortunately this problem can be fixed in a straightforward way. First we 
rewrite (12211 as 



#1 



1 



S< e iTX 



1 + A 



X I - 



Jo(Alog-) \ + ^{e 



A 



Ri 



du-J (A log 
it V A/ 



1) 



J, 



o (A log I) 



In order to study the last integral, we choose f(u) as above and h(u) = X[o,\]( u )i e ~ lTU ~ !]• The Mellin transform 
of /i can be written in terms of confluent hypergeometric functions but one can work with a simpler expression by 
realizing that, since f(u) is zero if u > A, one can replace h with a new function /i that differs from it only for u > A: 

-' irs - 1 s e (0, A) 

s e [A,oo) ' 



The result is easily seen to be independent of this extension of h. The integral giving the Mellin transform of h converges 



in the strip -1 < 5R(z) < 1 with M[h; z] = T^ z T(z)e 



-i-Kz/2 



A* 



Since we now have overlapping convergence strips 



we can use the Mellin-Parseval formula and obtain a Mellin-Barnes representation of our integral as 

t \- z T(z)e- l7TZ / 2 ~ i' 



1 

Ry 



^{ e iTX 



A 
2~7ri 



dz- 



VA 2 + z 1 



Ri 





u 


--[ 






2TTiJ c 



c+100 dz T- z \- z T(z)e-^ z / 2 



with 5R(c) e (—1, 0). We have employed the fact that 



dz- 



dz 



Va 2 



VA 2 + z 2 



0. 



(23) 



The analytic continuation of the integrand in the r.h.s. of Ij23(l is immediate. It has algebraic singularities at z = ±i\ 
(see fig. |SJ) with cuts in those points of the imaginary axis where \z\ > A, and simple poles on — z G N U {0}. 
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p = 



6 



/A 



r. 



7i 



72 



-i *A(p+l) 



73 



75 



76 



77 



iA|p-l| 



-iA|p-l| 



-o -iA(p + l) 



FIG. 6: Cuts and singularities of the integrand and integration contours in the complex plane for the Mellin-Barnes represen- 
tation used in the asymptotic analysis for r — > oo in the cases p = and p 7^ 0. 



In order to study the limit r — * + we can displace the integration contour in (|23|l leftwards parallel to the imaginary 
axis 8 . A simple analysis using the asymptotics of T(z) for large values of $s(z) shows that the series given by the 
residues of the integrand at the poles — z € N, namely 

Ri I ^fc!VA 2 + fc 2 y 

converges to the value of the integral. This is precisely what one would get by expanding the exponential in i|22|l . 
exchanging sum and integration, and computing the integrals that appear. This was the procedure used in [fl- 
it is not possible to directly get the asymptotic behavior in the r — > 00 limit by displacing the integration contour 
rightwards because of the presence of the cut and, even worse, the absence of poles with 5R(z) > 0. One could consider 
choosing a value for c and write (|23[) as a real integral in the variable y 6 K after the change z = c + iy. In fact, by 
doing that one obtains the sum of two ^.-transforms with asymptotic parameter logr that can be studied with the 
standard formulas. Unfortunately, proceeding in this way one gets a trivial zero asymptotic expansion. The reason, 
as we will find out later on, is that logr is almost, but not quite, the right asymptotic parameter 9 . 

The solution to our problem is displacing the integration contour all the way to the cut (choosing c — > 0~ ) as shown 
in fig. HJ The contribution of the arcs T2 and to the integral in the r.h.s. of i12.'SI) goes to zero as — c = e — > + . 
The remaining three contributions to the integral (without prefactors) are 

r ' ~j x * • 



Formula 1201 is obtained, precisely, by displacing the integration contour parallel to the imaginary axis. For functions with the asymptotic 

behaviors considered in the theorem the only singularities are poles whose residues give the asymptotic expansion. 

More precisely, an asymptotic sequence given by inverse powers of log t is appropriate to capture the behavior of our integral in r. 
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TTl „ 

- — +^P.V. 
A 



A e -iylog(Ar) e7 r !/ /2 r / i X 



— A 



-A e -iy log(Ar) e 7rj//2 r ^ y ) 



!J 



A 2 



The asymptotic expansions of these integrals can be obtained by a variety of methods, for example by steepest 
descents. Nonetheless, it is now possible to employ formula l|20(l and obtain, in principle, the asymptotics to any 
order. 

The integral over I\ can be written as the following standard ^.-transform by the change of variables y = A(l + s) 
with s £ [0, oo): 



_ e -iAlog(Ar) e 7rA/2 



ds e 



-iXs lot 



. T e-* As lo s A e 7rAs / 2 r[a(i + s)] 

y/s(s + 2) ' 



Taking = e~ 4Az , /(z) = 



iXz log A iii/2j 



— r[»A(i+z)] ^ an( j aS y m pt ^j c parameter logr we can use formulas (1191) 

y/z(z + 2) 

and H2()[). Since the asymptotic behavior of / when z — > + is 

/(z) ^£M + 0(z i/ 2) 



v2z 

and M[h\ z] = \~ z e ~T(z), we get for the integral over Ti an asymptotic expansion in powers of 1/%/logT with a r 
oscillating factor e - l Aiog(A-r)^ name iy 



_. iA log(Ar) + ^-^ r( iX) f /Z^Z + 

[ V 2Alogr 



O 



1 



(logr) 3 / 2 

A similar argument can be used to calculate the asymptotic expansion of the integral over T§, 

1 



JX log(Ar)- 



— + -T(-iX) 



2Alogr 



O 



(logr) 3 / 2 



Finally, to get the asymptotics of the integral over I3 we express it in the form 



+ i 



dy cos[y log : 



dysin[y log r] 



7^7 



By introducing neutralizers as in [l4] | it is possible to show that the critical points for the previous two integrals are 
and A. The contributions to the asymptotic expansion from both points can be obtained in a straightforward way 
using formula (|20|l . The contribution from zero is — ^ whereas the one from A coincides precisely with the sum of 
the contributions obtained before for Ti and IV Summing all the terms and substituting the result in wc finally 
get the behavior in r presented in Section ITD 

It is worth remarking that the asymptotic expansions discussed can be written as an oscillating factor that depends 
on logr multiplying an expansion in terms of inverse powers of y'logr. Moreover, note that the frequency of these 
oscillations is exactly that of the oscillating factor that we have already obtained. 



2. p^O 

Consider now the integral JSTJ). We set h(u) = e~ lru and /(«) = X[a,x]( u )i J o( Alo S j) •M/oAlog f). The Mellin 
transform of / is 



M[f;z] 



y.-2 
( 



A 2 (l + p 2 ) + (z-l) 2 
2pA 2 
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Proceeding as in the p 
1 



case we get the following Mellin-Barnes representation: 



1 



irX 



Q-i 



l 



-.s< e 



irX 



Q-i 



n 



i 

2nd 



2p 



1 

2nd 



c+ioo 

dz 

c—ioc 

c+ioo 

dz T 



iirz/2 

Z 



T- z \- z T(z)e 

" z / 2 T{z)Q_ 



Q-l 
A 2 (l 



A 2 (l + P 2 ) + z ; 
2pX 2 



o 2 )+z 2 



2p\ 2 



(24) 



The main difference with the previous case is the appearance of the associated Legendre function Q_ i . This function 
has cuts on the imaginary axis for > MP ~ M- O n the left of the cut the real and imaginary parts behave as 

shown in fig. {7\ They become singular in the neighborhood of the four points z = ±iX\p — 1| and z — ±iX(p + 1). 
The real functions that describe each of the continuous pieces are shown in the plots. Some of them are given by 
Q-i (x) evaluated on x > 1. Quite surprisingly, the others can be written in terms of the associated Legendre function 

P_i [P_i(x) — F (i, i; 1; ^5^)] with argument x > —1. The fundamental consequence of this is a change in the 
singularity structure of the integrand of l|24|l . The limit r — > + can be obtained as we did for p = 0. The result is 
given by (|TT)|l . 

By displacing the integration contour rightwards to the cut as in the p — case (see fig. 01 we can split the last 
integral in l|24|) into seven different contributions corresponding to the integration curves 71, ... ,77. It is straightfor- 
ward to show that the contributions from 73 and 75 go to zero as the contour gets closer to the imaginary axis. We 
are then left with 



7i 



72 



74 



76 



77 



-iy log(A-r; 



e*y/ 2 T(iy)Q_i 



dy e 

A(P+1) 
A(P+1) 

dy e-^^e^Tiiy) 

X\p-1\ 

2 



P_ 



1 



2p 



iP.V. 



X\p-1\ 
dy e 

-X\p-1\ 



y 2 -X 2 (l + P 2 ) 
2pX 2 

y 2 -X 2 {l + P 2 ) 
2pX 2 

iylo S (\T) e *y/2 T ( iy jQ_ 



iP_ 



7T 
— I 

2 



-a|p-i| r 
dy e-^ log(AT) e Tra/2 r(i?/) P_ 

X(p+1) 
-A(p+1) 

dy e-^ los(AT) e Tra/2 r(i?/)Q_, 



y 2 -X 2 (l+p 2 ) 
2pX 2 

< 2 -X 2 {l + p 2 ) 
2 P X 2 



-iP_ 



X 2 (l + p 2 )-y 2 
2pX 2 

X 2 (l + p 2 )-y 2 
2pX 2 

X 2 (l + p 2 )-y 2 
2 P X 2 



All these integrals can be written as Fourier transforms whose only critical points are the integration limits. These 
can be isolated by using the appropriate neutralizers ^4[ and the integrals so obtained can be written as /i-transforms 
by straightforward changes of variables. An interesting feature of some of the integrals that arise is the fact that the 
asymptotic expansion of the corresponding / functions involve logarithms precisely in the form assumed in Q18f) . This 
means that in addition to the p m o coefficients that appeared in the p = case we will have also contributions coming 
from the p mn 's with n^O. 
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In order to get the asymptotic behavior of the different / functions close to their singularities it is useful to remember 
the asymptotic expansions of P_i(x) and Q_i (x) at x = — 1 + and x — 1 + , respectively: 

5 1 

P_ i (x) ~ - log 2 log(x + 1) H 

2 7T 7T 



Q_i(ar)~-log2--log(^-l) + --- 

With the previous guidelines, and following the same steps as in the p = case, it is easy (albeit lengthy) to 
deduce @. It is possible to arrive at the same result by different methods that are somewhat simpler (i.e. steepest 
descents) but they only allow to compute the first contribution to the asymptotic behavior of the field commutator. 
The advantage of the method employed here is that one can calculate as many terms as desired for the asymptotic 
expansion just by taking more terms in the asymptotic series of the function /. 



C. Asymptotic behavior in A 



We will consider now integrals of the form 

i-A 



B 



dq f dh f A 2 / { o)(g;ii,t 2 )e A * (a5tl ' t » ) 

'7i ^72 



(25) 



(0 < B < A) where q is a real variable, t-y and t 2 are complex variables and 7i )2 are closed paths (possibly different) 
in the complex plane. The functions /( ) and $ are supposed to be holomorphic in the two complex arguments and 
C n in q (with n possibly infinite) in an open neighborhood of the integration region. In the following we will use the 
notation d { = d u , i = 1, 2, V = (d q ,d 1 ,d 2 ), and ||V$|| 2 = i] q (d q ^) 2 + i] tl (d 1 ^) 2 + 7ft 2 (d 2 $) 2 , where r) q , r) tl , and r/ t2 
are some functions of q, t%, and t 2 that will be chosen later in the discussion. This freedom will prove to be useful 
when we parametrize the curves 7i i2 because it will allow us to have denominators in the integrand that are real 
positive functions, facilitating the analysis of the A — * oo limit. If || V<I>|| 2 differs from zero for all q 6 [B, A], t\ £ 71, 
and t 2 € 72, it is possible to devise an "integration by parts" procedure for this type of integrals by writing them as 



dq <s> dt\ (t dt% 
J 71 ^72 



f i0) v q d q ^d q e^ / (0 )^i9i1>9ie A * f^Vt^^e 



|V$| 



The contribution of the first term can be written, after integrating by parts in q, as 

A 



dt\ i> dt-2 
71 ^72 



f(o)Vqd q ^e 



A4> 



B 



dt\ i> dt 2 

71 ^72 



A 



dqe^d, 



B 



f(o)V q d q <S> 
' IIV$|I 2 



(notice that the integration region is compact and the function integrable so that we can apply Fubini's theorem and 
freely exchange orders of integration). Integration by parts in t\ gives 



dq d> dt 2 f dt\d\ 

-'72 71 



/(O)^ti^i^ e 



dq (p dt 2 f dtie^di 

-'72 •'71 



/(O)^i^l^ 

liv$n 2 ■ 



If (_f(o)f?ti<9i<I>e A *)|| V^H -2 is a holomorphic function of t\ in an open set containing the closed curve 71 (for all 
q G [B, A] and t 2 G 72) the first integral is zero and we are left only with the second term. The same argument applies 
to the integration in t 2 . We finally get that l|25|l can be written as 



-xi dti 



dt 



72 



dq f dh f dt 2 f {1) (q-t 1 M)e mq ' MM \ 
'71 J 12 



(26) 



here we have introduced the notation 

_ „ f(k)V q d q t> 



/(fc+l) = d q - 



|V$|I 2 



f(k)Vtidi® , g f(k)Vt 2 d 2 <f> 



|V$| 



|V$| 



fc = 0, 1,. 
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The second integral in <|26(l has the same structure as the original integral, so the same procedure can be applied 
as long as the appropriate regularity conditions are satisfied for expressions involving the successive k — 1, .. . 
The strategy to obtain asymptotic expansions for the original integral consists in getting expansions for the lower 
dimensional boundary integrals that arise in this way, provided that the remaining non-boundary integral behaves 
nicely in the A — > oo limit. Particularizing this argument to the case of a single complex variable t is straightforward. 

In our case [see equation Hl()|l ] we actually want to obtain an asymptotic expansion of an improper integral. So we 
will also have to study the asymptotics of the limit A — > oo, 




In order to derive useful asymptotic expansions from the previous formulas it is desirable that e A *( 9; * 1 '* 2 ^ be bounded 
by a A-independent function 10 . This happens for example if !ft[<I>(q; ti, £2)] < in the integration region [B, A] x 71 x j 2 - 
In fact, we will take advantage of the freedom in the choice of the integration contours 7x,2 to impose this condition. 

If there are points in the integration region (referred to in the following as singular points) where ]|y$|| 2 = the 
argument presented above to derive (j26|) is no longer true. In such a situation the way to proceed is to isolate 

them by introducing neutralizers 11 . This allows us to divide the integration region into several pieces and study the 
resulting integrals by choosing the most suitable techniques. This is facilitated by the fact that neutralizers force most 
of the boundary integrals to be zero. In the process of obtaining asymptotic expansions from Ij27(l it may be useful 
to further divide the integration region of the last integral in several pieces, by introducing additional neutralizers, in 
order to localize the singular points in a region that does not contain the boundary q = A; this simplifies the analysis 
of the A — > 00 limit. 

The freedom to choose the integration contours in (|25|) may be useful for other purposes. One can, for example, 
avoid singularities of the boundary integrals over {£>} x 71 x 72 and {A} x 71 x 72; it may be possible to write the 
integrals involving the singular points as (multidimensional) Laplace or Fourier transforms for which the analysis of 
singular points can be carried out by standard procedures (see 14] ) or to improve the behavior for large values of q 
if one intends to study the asymptotics of the improper integral when A — > oo 12 . Finally in some cases this freedom 
may allow us to simplify the analysis of the contribution of some of the singular points. 

As in the study of the r asymptotics, we will consider the cases p = and p 7^ separately. 



1. p = 

We study the integral lfTU)l for p = 0: 

Let us discuss first the possible choices of integration contour 7. As commented above it is convenient to impose that 
the real part of the exponent in the integrand be less or equal to zero, i.e. — 1/t) < 0. It is straightforward to 
show that this happens when u(u 2 + v 2 — 1) > (t = u + iv). This complex region (see fig. |BJ is given by the points 
inside the unit circle centered at the origin with positive real part, those outside this circle with negative real part, 
and the boundary of these two regions (i.e. the imaginary axis and the unit circumference centered at the origin). 
The integration contour can be any closed, positively oriented, simple curve contained in this region that surrounds 
the origin. Notice that all these contours contain the points t = ±i. If we choose rj q = — 1 and rj t = t 2 we have 



The asymptotic analysis of multiple integrals is greatly simplified by the identification of the critical points, those points that give the 
dominant contributions. For Laplace or Fourier types of integrals these are easy to identify and, in practice, they are just a finite number 
of isolated points that can be singled out and studied by using neutralizers. 

These are real, positive, C7°°(M) functions v(q\ ct\, 02; /3i, $2) satisfying v(q\ ct\,a2\ P\, P2) = if q £ (— 00, ai] U [02,00) and 
v(q\ «i , ol2\ /3i, (82) = 1 if q £ P2] (some of the parameters oji , 02, /9i, 02 may be taken to be infinite). In several variables q\, . . .q n 
it is sometimes useful to consider neutralizers that depend only on \\q\\. 

It is important to point out that the asymptotic behavior of an improper integral over [0, 00) may be very different from the limit of 
the asymptotic expansion of the integral over [0, A] when A —* 00. 
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|V$| 



t + 



- It- 



2ire 



At the left boundary (q = 0) this is zero when t = — ir±v 1 — t 2 . Depending on the value of r these points may be on 
the unit circumference; it is, hence, convenient to choose integration contours in the allowed part of C that avoid them 
(see fig. |SJ|. This will always be possible except for r = 1. When q ^ the condition ||V$|| 2 = gives, generically, 
four possible solutions (depending on q and r) that may or may not be in the integration region. If r > 1, there is one 
of these, (t = —i, q = logr) that will be in the integration region for every possible choice of the integration contour 7. 
The remaining ones will or will not belong to the integration region depending on 7. In principle one would proceed 
by picking a certain contour, determine the values of q corresponding to these points, introduce suitable neutralizers, 
and study the neutralized integrals one at a time. There is, however, a simple argument that shows that it is not 
necessary to discuss them in detail (we assume q 7^ in the following discussion unless otherwise specified). 

The idea is to write the integral in (|28|l as the sum of three integrals Ij, j — 1, 2, 3 and choose different 13 contours 
for each of them. These integrals are 

where we have introduced three neutralizer functions Vj{q) : j = 1, 2, 3 satisfying v\ + V2 + ^3 = 1 in [0, A] and 

v\{q) = 1 if q e [0, a>i], 



The reason why we use neutralizers depending only on q is, precisely, to exploit this freedom. 
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= 


if 


q e ["2, A\, 


V2{q) 


= 


if 


q e [o,ai] u [/?2,i4] 


V2(q) 


= 1 


if 


q e [a2, A], 




= 


if 


?€ [0,/3i], 


^3(9) 


= 1 


if 





with < ai < a2 < /3i < P2 < A (the specific choices for these parameters will be explained later). By doing this the 
effective integration regions in q are [0,0:2], [oi,/?2] and \J3±, A] respectively; notice that the boundary q = appears 
only in the first. 

Let us consider first I\. In this case it is convenient to choose an integration contour 7 that avoids the unit 
circumference (except at t — ±i) because, if t < 1, the integration by parts procedure introduced above could, 
otherwise, give a boundary contribution at q — involving a singular function. We pick the contour depicted in fig. 
[SJ In order to use integration by parts it is also necessary to show that it is possible to find 02 such that ||V<1>|| 2 7^ 
in the integration region. The solutions to ||V$|| 2 = are those satisfying the following quadratic equations 

(1 + q)t 2 + 2iTe~ q t - 1 + q = 0, 
(1 - q)t 2 + 2iTe~H - 1 - q = 0. 

It can easily be shown (by using the implicit function theorem or the explicit solution of these equations in t) that 
the solutions are continuous at q = (if r 7^ 1), i.e. the solutions for small enough q are close to t = —it ± yl — t 2 . 
If r 7^ 1 it is then possible to choose 02 and 7 in such a way that || V^lj 2 7^ in the integration region 7 x [0, 02]. 

We can use formula (|27[> to get an asymptotic expansion for 1\ . We first note that the term corresponding to the 
second integral in l|2t)|) is 0(1/ X) because our choice of contour allows us to bound the absolute value of the integrand 
by a regular function independent of A. The integral I\ is then given at leading order by the contribution of the 
boundary terms which can be written as 



h 



dt 



7ri?i 1 Lt 2 + 2irt - 1 



[the function v\(q) kills the contribution coming from the q = 02 boundary]. If r < 1 the poles of the integrand 
t = -IT ± Vl - T 2 are on the unit circumference but only that with the negative real part is inside 7. The value of 
the residue of the integrand at this point is — 2v / 1 1 _ T a and hence 



If r > 1 the poles of the integrand are on the imaginary axis t = —it ± «V r 2 — 1; now only that with the positive 
sign in the square root is inside 7. The residue is — A , and hence 

h = wA7^t] = l^f (29) 

Consider next 12- In this case it is convenient to choose the unit circumference for 7 because, parameterizing it by an 
angle (t(9) — e 10 , 6 € (— tt, tt]), the exponent A$(g; t) becomes a purely imaginary function of the two real variables q 
and 9, 

iX^O) = \$[q;t(6)] = iX(qsin9 - re" 9 ), 

and the integral becomes a Fourier transform on a rectangular region determined by the support of the neutralizer. 
As is well known the critical points for Fourier transforms are given by the boundaries of the integration region and 
those interior points where the gradient of the real function <f> is zero. In our case the q — a% and q — P2 boundaries 
give no contribution for this integral owing to the presence of the neutralizer. Let us then find the interior critical 
points. Since || V$|| 2 [g; t{9)] is a sum of squares of real numbers 1 these points are precisely those where our procedure 



This is no longer true for other choices of 7 for which || V$|| 2 is, in general, complex. In this cases 1 1 X7<E> 1 1 2 may be zero even though the 
partial derivatives of $ are different from zero. 
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of integration by parts fails. We have to solve the following equations 

qcos9 = 0, 



■ re 



0. 



The solution is q = 0, 9 = — arcsinr if r < 1 (these are outside the effective integration region for I 2 ) and q — logr, 
9 = — ^ if r > 1. The choices of can be made in such a way that the critical point with q > (present only when 
t > 1) always corresponds to I 2 (we also choose /3i such that logr < 0x to avoid having critical points in The 
asymptotic expansion of 1 2 can be now obtained by using the standard formulas for the contribution of critical 
points. In our case this is simply 



Ri 



= iA(t — logr — 1) 



VlogT 





1" 






A 



(30) 



for r > 1 and if r < 1. We do not need to discuss the limit A — > 00 of I\ and I 2 because the integrand has compact 
support independent of A; we have only assumed that log r < A. However this is no longer true for I 3 . 

Let us analyze, finally, lim^^oo 23 . We choose again the unit circumference for 7 (with the same parametrization) 
and write 



lim .Z3 

A^oo 



Xe irX ,. 
- — — lim 

27ri?i A— >00 J a 



dq / dS iy 3 (q)e tX ^ q -' 9) 



Using formula (|27|l with r\ q 
Xe lrX 



T, r\t = t 2 , and taking into account that ^(Pi) = 0, this expression becomes 



lim 



X 



d6 



(re~ A + 8^0)6 



i\<t>{A;e) 



A 2 , 



' + (re- A + 1 



dq d9v 3 (qy X ^ gi (q;9) 



(31) 



Here 
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gi {q-9) = ie M f (1) [q;t 



(rer q + s\n9) 2 (Ter q - gsinfl) + q cos 2 9[(q 2 - 4) sinfl - r(q + 4)e~ 9 ] 
[q 2 cos 2 9 + (re-9 +sin6>) 2 ] 2 



We have disregarded in equation (|31|l a term containing derivatives of v%{q) because it can be shown to be 0(1/ X R ) for 
every R > 0. We want to show that (|3*T|l does not contribute to the asymptotic expansion at the order in A considered. 
Note that it is difficult to prove that the limit of the last integral is zero when A — ^ 00 by using an argument inspired 
in the Riemann-Lebesgue lemma, because changing variables to write it as a Fourier transform introduces singularities 
in the integrand that are not easy to deal with. To study the asymptotics of (|31[) we will follow instead several steps: 

i) Use integration by parts to further decompose the second integral in (|31|l as a surface term and a double integral 
in q and 9. 

ii) Split each of the integrands obtained in this way in two pieces; one with a simpler denominator and another whose 
limit A — > 00 vanishes. After this step we will be left with two integrals in q and a double integral. 

iii) Split in two the double integral by writing the exponential term as 

iv) Show that the contribution coming from e lA 9 sme ( e -* A ' re ' — 1) vanishes sufficiently fast when A — > oo by using 
the Holder inequality. 

v) Finally, show that the remaining terms have a simple asymptotic behavior. 
Step i) After an additional integration by parts l|31ll becomes 



Xe iTX ,. 
- — — lim 

27riti A^oo 



X 

-i\re~ 



d9- 

n A 2 cos 2 9 + (re- A + sin 9) 2 
{sm9-Te- A )e lXAsine 



d9- 



X 2 7_ 7r A 2 cos 2 9 + (re~ A + sin 9) 2 
J_ [ A 

A 2 4 



9i(A;8) 



dq / d9v 3 ( q y x ^g 2 (q;9) 



15 In the following we will denote gk(l', S) = *e l8 /(fe) the «e l6 factor comes from the measure in the integral. 
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Step ii) All the integrands that appear in the first step have factors of the type 

1 



[g 2 cos 2 9+ (re-i + sin6>) 2 ] 



N 



N = l,2, 



that can be written as 



Z)j=i T 3 e^ jq Pj(q, sm9, cos 9) 



+ 



[q 2 cos 2 9 + sin 2 9} N ' [q 2 cos 2 9 + {re-i +sm9) 2 ] N [q 2 cos 2 9 + sin 2 9} N 



N = 1, 2,... 



with Pj(g, sin#, cos 0) a polynomial of degrees 2 (TV — 1), 2 (TV — 1), and 2N — 1 in q, cos 9, and sin# respectively. For 
fixed values of q larger than a certain Q > logr the maximum of the function y{9) = [q 2 cos 2 9 + (re~ g + sin^) 2 ] -1 is 
(1 - Te~ q )~ 2 , so we have 16 



Ei=i T j e-™Pj (q, sin 6, cos 0) 



[q 2 cos 2 9 + {re-i + sinfl) 2 ]^ 2 cos 2 9 + sin 



n 2 fflN 



< 



(1 -re~ q ) 2N {q 2 cos 2 9 + si: 



, 2 f)W 



N = 1, 2,... 



where -P,(g) is the polynomial in g obtained from Pj(g, sin#, cos#) by taking the absolute value of the coefficients and 
substituting sin6> = 1, cos# = 1. This formula, and the fact that 



d6- 



2 4 r(i !-. 1/2) asA^oo, 



_ w [A 2 cos 2 + sin 2 9] N ' " A T{N) 
allows us to show that the first two integrals obtained in step i) can be substituted by 

sin# 



d9 



e >XAsine + ± r d6 MA 2 - 4) sin 2 9 cos 2 9 - A sin 4 9 ^ A sin e ^ 



A 2 cos 2 9 + sin 2 9 ' A 2 J_„ [A 2 cos 2 9 + sin 2 9} 3 



Let us consider now the double integral appearing in step i). The function 52(9; 0) can be written as the sum of two 
pieces h 2 {q;9) and G 2 (q;9). The first one is 

h 2 (q;0) = 

g 4 (g 2 - 4) cos 6 9 + g 2 (40 - 23g 2 + 3g 4 ) cos 4 9 sin 2 9 - (4 - 23g 2 + 8g 4 ) cos 2 9 sin 4 + (g 2 - 1) sin 6 
* (g 2 cos 2 6> + sin 2 (9) 4 

and G 2 (g; 0) is a sum of factors {re~ q ) N (N = 1, 2, . . .) times quotients of polynomials in g, sin6>, and cos9 divided 
by products of powers of [g 2 cos 2 9 + sin 2 9] and [g 2 cos 2 9 + (re~ q sin 9) 2 ] (notice that these functions never vanish in 
the integration region). The integral becomes 

-72 t d( i r d9v 3 (q)e iX ^[h 2 (q;6)+G 2 (q;9)]. 

A J 01 J-TT 

The absolute value of the integrand involving G 2 (q;9) is an integrable function owing to the decreasing exponential 
factors and the fact that they multiply terms that grow polynomially at worst in g. Hence its contribution to the 
total integral is 0(1/A 2 ) even after taking the limit A — > 00. We thus conclude that 



Xe 



lim 
27ri?i A^oo 



sin# 



i\A sin 6 



-if- 

A J_ n A 2 cos 2 9 + sin 2 9 
1_ r dQ A{A 2 -4) sin 2 0( 
+ A 2 J_ v [A 2 cos 2 9 + sin 2 9] 3 

~~^C dq f deelUiq '' 9)v ^) h2 ^ 6 ) 



de A{A 2 - 4) sin 2 9 cos 2 9 - A sin 4 9 &iXA sin 9 



(32) 





1" 






A 



16 We choose (5\ > Q > log t in the neutralizers introduced before. 
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Step iii) Note that we have already managed to substitute the exponential factor e lX <t>( A 'fi) by e lAsine within the 
integrals over 9. However it is not possible to remove the function <f>(q, 9) from the double integral. We instead 
proceed to rewrite it as 

-T2 t d 1 I' d6e iX * sine [e- iXTe -* - l]v 3 (q)h 2 (q; 6) - ± [* d q f d9e* Xqsine v 3 (q)h 2 (q;9). 

A J0! J-TT A J/3 1 J—rr 



Step iv) Consider the first of these integrals in the limit A — > oo, 



1 

A 2 



dq / d9e lXqsin6 [e- lXTe ~ q - l]u 3 (q)h 2 (q;9). 



(33) 



This integral converges because (e lXre q — 1) falls off exponentially as q — * oo. Let us define now the functions 
f(q,9) = q°{er iXTerq - 1) and g(q,9) = q- a e iXqs ™ B i> 3 {q)h 2 {q]9) with a chosen so that f,g e L 2 ([^,oo) x [-tt.tt)). 
For example we set a = 3. The Holder inequality gives 



1 

A 2 

< 



A 2 



dq I d9e lXqs ' me [e- lXTe 
dq f d9\f(q,9)\ 2 



l]u 3 (q)h 2 (q;e) 



< 



dq d0\f(q,0)g(q,l 



dq / d9\g(q,l 



The integral involving \g(q, 9)\ 2 is convergent and independent of A and the one involving \ f(q,9)\ 2 satisfies 



dq I dd\f(q,9)\ 2 y <2V2^\ I dqq e sin 2 ( ye"« 



2V2tt<{ / fft-(log£) 6 sin 2 (Ai) 



after changing variables (t = e q ) and defining A = tA/2. As we show in Appendix I this last integral is 0(log 7 A) 
and, hence, the term (O is 0(log 7/2 A/A 2 ). So in expression I|.'i2[| we obtain 



_5 



Ae iTA ,. 
hm 

2ttRi A— »oo 



' A 
1 

A 2 



d6»- 



^■iAyl sin ^ 



sin 2 6» 



ig ^4(A 2 - 4) sin 2 9 cos 2 g - A sin 4 g sin e 



[A 2 cos 2 6» + sin 2 6>] 3 



(34) 



— 7T 

,4 





"1" 




riog 7/2 Ai 


} + o 




+ 


A 


A 



The reason why we had to integrate by parts in the first step is to get an additional factor of A dividing the powers 
of log A. 

Step v) Using integration by parts twice one can check that the terms in the square brackets in Q34JI can be written 
when A — > oo as 



2tt lim 



dqu 3 (q)J (Xq)= lim / dq d9 is 3 (q)e lXqsine , 
A— too 



(35) 



plus an extra contribution coming from derivatives of v 3 (q) that does not contribute at this asymptotic order. In 
Appendix II we show that is 0{l/\ R ) for all R > 0. We therefore conclude that (JSTJ) is 0(log 7/2 A/A). In our 
proof we have used integration by parts twice to arrive at 134(1 . Actually by using it repeatedly and applying the 
five-steps procedure explained above it is possible to argue that limA^oo ^3 is, in fact, 0(1/ X R ) for all R > 0. 

We hence conclude that the first terms in the asymptotic expansion of l|28|l are given by the contributions 129fl and 
(|30|l . and therefore by when r > 1. In order to get the first non-zero term in the asymptotics for r < 1 it is 
necessary to perform integration by parts twice and follow the steps detailed above. This leads to the contribution 
for t < 1 shown in l|12fl . 

1 (i f Strt 2 \ = 1 r(l + 2r 2 ) 

2ttR 1 \ A J {t 2 + 2irt - l) 3 J R x 2A(1 - r 2 )5/ 2 ' 
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The A — > + limit, on the other hand, can be obtained by using Taylor's theorem to write 



E 

fc=0 



(N + iy. 



with < £(t) < r\e~ t/x . Substituting this in © (with p = 0) we get 



1 



e iTX 



dtJ (t) 



N 

E 

k=Q 



(~irX) k _kt | (— irA) 



(JV+l) 



fc! 



(JV + l)! 



The contribution of the last term can be easily bounded 



e' irA HrA) 



JV+l roo 



(jv + i)! y 



< 



(rA) 



JV+l /-OQ 



(JV+1)!7 



r^A^ 2 

cite — r-* = 



(JV + 1)(JV + 1)! 



= 0(X N + 2 ) 



So we find that the asymptotics of in this limit is given by 

JV 



A. 

Ri 



rAV . (-irX) k 
kWk 2 + A 2 



0(A 



JV+2-, 



We study now the integral 1)1 0|l for p 7^ 0. We will basically follow the same steps of the p — case so we will skip 
some details. We choose rj q = — 1, r) tl = tf, and = i 2 , obtaining 



ti + 



p 2 1*2 



ti - — +p(t 2 - — ) + 2i? 
*1 £2 



As we did above we introduce neutralizers ia,-, 7 = 1, 2, 3 and write (|10|) as a sum of the tree integrals 



h = - 



Ae 



4tt 2 



•I 71 J ~I2 



ht 2 



-)— ire 



Starting with ii we fix 7^2 as in fig. [S] Integrating by parts we find that only the boundary term contributes at 
leading order. This contribution can be written as 



* 2 



1 



, :j pt 1 (i 2 -l)+t 2 (i 2 + 2zrt 1 -l) 



(36) 



This integral can be computed exactly (see Appendix III) and, in fact, it is equal to the free commutator given in J6(]; 
in particular it is zero in the regions IA and IB of fig. This means that we will have to use integration by parts 
again, as in the p = case, to get the first significant asymptotic term when A — > 00 for these regions. After doing 
this we get the double integral 



tit 



1 '2 



., .,- 2 [pt l (t 2 2 -l)+t 2 (t 2 1 +2iTt 1 -l)} 3 



(37) 



It can be easily checked that this reproduces the result obtained above for p = 0. This integral, in the considered 
regions outside the light cone, can be written in terms of complete elliptic integrals of the first and second kind, as 
shown in Appendix IV. The result is just the 1/A contribution appearing in equation (|13J) . 

Though it is also possible to compute the integral I\ in the remaining regions, we will not do so because their 
contributions are subdominant with respect to those coming from the critical points of I 2 . 

Consider then I 2 and take as before the unit circumference as the integration path 712- We parameterize this curve 
according to ti{0\) = e lSl , t 2 (9 2 ) = e 1 ® 2 , 1)2 G (-7r,7r]. Now we have 



i\<j>(q;l 



2) = A$[g;*i(0i),i 2 (02)] = iA[g(sin0i + psm0 2 ) - re~ 
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so that the critical points are given by the solutions to the equations 

sin Oi + p sin 62 + re~ q = 0, 
q cos 6*i = 0, 
qp cos 82 = 0. 



Since we are computing I2 we are only interested in solutions with q 7^ 0. Therefore we must have cos 6\ — and 
cos 62 — 0, i.e. 61 = ±^ and 82 = ±§. For 6\ = 62 = ■§, the remaining equation implies 1 + p + re~ q = 0, which 
has no solutions for 9 e I. If 8\ = —62 = § we get 1 — p + TeT q = 0; this equation has solutions q — log -33- in 
the integration region q > only if r > p — 1 > 0. If 8\ = —62 = — f-, we must have — 1 + p + re~ q = 0, and there 
exists a critical point q = log in the integration region only if r > 1 — p > 0. Finally if Q\ = 62 = — \ we obtain 
— 1 — /? + re~ 9 = 0, which is solved by q = log j^— ; in this case g > only when r > 1 + p. As we see when (r, p) are in 
the regions IA and IB of fig. ^there are no critical points for I2, if (r, p) is in region II there is only one critical point, 
and if (r, p) is in region III there are two critical points. The contribution of these critical points to the asymptotics 
of I2, that can be obtained with the standard formulas [T^J, are the following in regions II and III respectively: 



-a 



-if^MT+lp-lKl + logT-^r)] 



Ri [ y/2n\p\l - p\ log 





1 


}+0 


_A 3 /2_ 



Ri 



-if p iA[T-|p-l|(l+log j^jj)] gi f el A[T + (p+l)(log i±*-l)] 



V2tAH1-p| log ^ V2tV(1 + p) log 





1 


} + o 


_A 3 /2_ 



These, together with the contribution of the boundary term, provide l|14[) and l|15[l. 

Finally it is possible to prove that lim^^oo ^3 gives no contribution by essentially following the same steps as when 
P = 0. 

It is interesting to comment on the singularities that appear in the borders between the different regions I, II, III, 
and the axis p = 0. Their presence is a manifestation of the fact that the asymptotic behavior in A changes abruptly 
between adjacent regions. Notice, nonetheless, that the behavior is continuous in the border between the portion of 
the axis with r < 1 and region IA, where the leading terms of the asymptotic expansion in A are both 0(1/A). 

To conclude let us point out that the limit A — > + can be obtained as in the p — case. The result is 



1 



■KRi^/p 



r " (-»tA)\ 

fc=0 



A 2 (l + P 2 ) + fc 2 
2pA 2 



+ 0(X N+2 ). 



IV. CONCLUSIONS AND COMMENTS 



We have analyzed in detail the issue of microcausality for linearly polarized cylindrical waves by looking at the 
commutator of the axially symmetric scalar field that encodes the physical degrees of freedom of the system. We have 
been able to show several interesting effects that appear in the model. The first is a smearing of the cylindrical light 
cones. This is especially obvious if one studies the behavior of the commutator (divided by 8iG) at two points with 
coordinates (ti,R±) and (^2,^2) when A (the quotient of R\ and the Planck length) is large. If R2 is not too close 
to the axis one gets, in this limit, the discontinuous, cylindrical light cone structure defined by the free commutator. 
In particular, outside the region defined by this light cone the commutator is zero and, hence, observables would 
commute as in ordinary perturbative QFT. If, instead, one looks at the behavior when R2 = one finds a peculiar 
behavior as A — > 00; now superimposed to the free contribution, the commutator shows a characteristic oscillating 
behavior when the variable r = (£2 — ii)/4G grows. The frequency of this oscillation is controlled by the value of A 
but the amplitude is independent of it, and turns out to decrease very slowly with r in such a way that one recovers 
the value given by the free commutator only for very large values of r. Nonetheless, there is a sense in which the free 
propagator is recovered if one averages over time intervals sufficiently larger than AG. 

In our study we have had to find the most efficient methods to obtain the relevant asymptotic behaviors. Though 
this has not required the introduction of completely novel techniques it has been necessary to extend to our situation 
the usual Mellin transform methods and those employed for the analysis of multiple integrals. We have needed also 
to take into account that the integrals that define the commutator are, actually, improper. 
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Several issues are worth discussing in more detail. One is the problem of introducing regulators in the field operators 
and compare the results with those derived here. Though it is clear that the approximation provided by the improper 
integrals considered in this work must be good in certain regimes, at some point the existence of a physical cut-off 
might manifest itself in the behavior of the field commutator, especially in the asymptotic regimes that we have 
explored. Another task that can be confronted with the techniques that we have developed here is the computation of 
other matrix elements for the commutator. It would also be interesting to consider other Green functions and some 
related objects, such as S matrix elements. We plan to do that in the future. 



V. APPENDIX I 



We want to investigate the asymptotic behavior of J dt^^-— sin 2 At, m € N, when A — > oo. This can be found 
by a straightforward use of the Mellin-Parseval formula l|19f) in order to get a Mcllin-Barnes representation for the 
integral. Defining h(t) = sin 2 t and f(t) = ( log *) we get the Mellin transforms 



M[h; z] = - 



r(z) cos 



77 "_' 

2 



M[f; l-z] = - 



2 1+z 

r(m + 1; 



ytn+l 



-2<K(z) <0, 
1 < ft(z), 



so that 



I 



, (logt) m . 2 , ml 
d V 6 1 sin 2 At = — 
t 2m 



C+ZOO 



r(z)cos^ 

dzA- z ^ -2- ce (-2,0). 



The integrand in this last expression can be analytically extended as a meromorphic function to C with poles in 
z = — 2ri, ri = 0, 1, . . .. By displacing the integration contour to the right of the pole at z = 0, we get 



1 , (logi) m . 2 . 
d V 5 ' sin 2 At 



A" 



, mlT(z) cos 



2l+z z m+ 



ml f Xa+lco T(z) cos 
dzA 2 

2iri 



7T 

2 



2l + ^ r m +l 



(38) 



(xq > 0) whenever the integral converges and provided we can neglect the contributions from the segments needed to 
close the integration contour at large values of Q(z). Writing z — xq + iy we have 



ml f x "+ loo i r(*)cos 
dzA z 

2-ki 



2i+2 r m +i 



< 



4tt(2A) :e o J_ 



dy 



T(x + iy)cosf (x +iy) 



(x + iy) 



m+l 



and using lim 



-i=|r(a; + iy)\e^y\l 2 \y\- x+1 /' 2 = 1 we find that 



r(a; + iy) cos |(x + iy) 



% ix -m-3/2 



(x + iy) m+1 

Hence, the integral in l|38(l is absolutely convergent if xq < m + \ and its contribution is O [j^j] ■ It is straightforward 
to show that the residue at the pole z — is an m + 1 degree polynomial in log A with the higher degree term given 
by ( 2~(m+i) (logA) m+1 . This proves that, for m E N, 

..Qogt) 



dt- 



t 



sin 2 At = 0[(logA)(" l+1 )] 



VI. APPENDIX II 



Let us analyze the asymptotic behavior of 2tt dqis3(q)Jo(\q) when A — * oo. After changing variables according 
to k = Xq this integral becomes 

2-k r°° / k\ 
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Taking into account that Jo(k) 



Id d 2 

k dk dk z 



Jo(fc), this can be written as 



2-k f°° ( k 

— dkv 3 [ - 



A 



1 d d 2 
k dk dk 2 



N 



Jo(k) 



with JVgN. It is straightforward to prove by induction that 

N 2N 



1 d d 2 

k dk dk 2 



Ea 
U2h 



k 2N ~j dki 



for certain coefhcients <Zj G Z. Since all the derivatives of fc2 jv-j v z (t) at fc = A/3i cancel, by integrating by parts j 
times and changing variables according to k = qX we obtain for our integral the following expression 



2N 



k 

k 2N-j P3 [J 



9tt 

■*>(*) = w £(-i 



2iV 



X 2N 



Using that Jq(x) < 1 for all i e 1, we hence get 



27T f°° ( k\ 27T f°° 

— dkv i [jj J o(fc) < 3^]vEN J f3 dq 



d,J 



dq} 



A dj 
d<1 dql 



^3(9) 



^3(9) 



2N-j 



Jo(Xq)- 



2N-j 



with all the integrals in the last sum being convergent. We conclude that, for all N £ N, 



2tt r°° 

A Aft 



dkv z ( - ) Jo(fc) = 



1 



A 2JV 



VII. APPENDIX III 



In this appendix we compute the integral H^6(l . It is possible to show that choosing 7^2 as in fig. [HI its denominator 
is always different from zero except for those exceptional values of p and t corresponding to the borders between 
regions I, II, and III. We can use Fubini's theorem and exchange orders of integration. To integrate in t\ we fix a 
value of t 2 on 7. The poles of the integrand are 



1 

2h 



p-pti- 2irt 2 ± Jit 2 + (-p + pt\ + 2iTt 2 ) 2 



One of them is always inside 7 and the other outside. For a fixed value of t 2 we find that t 1 is inside 7 when 



If 



> 



or in certain exceptional cases when di(t 1 ) = 0. Analogously is inside 7 when 



j-\/4t 2 + (-p + pt 2 + 2iTt 2 ) 2 



< 



or, again, in certain exceptional occasions with = (in this last case t 1 is outside 7). The residues at these 

poles are 



tf 



JW 2 + (-p + pt 2 + 2irt 2 ) 2 ' 
-1 

J4t 2 + (-p + pt 2 + 2iTt 2 ) 2 ' 
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and integrating in t\, i|36|) becomes 



ttRi 



dt 



sgn R[| ^4t 2 + (-p + pt 2 + 2lrtp 
v/4t 2 + (-p + pt 2 + 2irt) 2 



(39) 



Let us define A = — 1 + p — r, £? = 1 + p + t, C = 1 + p — r, and D = — 1 + p + r. It is not difficult to check that 
4t 2 + (-p + pi 2 + 2irtf = -[A(t - if + B(t + if][C{t - if + D(t + if], 



and, hence, the integral 139JI is 

2i 



dt 



ignft[±y/[A(t - if + B(t + Q 2 ][C(£ - ^) 2 + D{t + i) 2 ] 
- *) 2 + S(t + if][C(t - Z ) 2 + £>(i + if] 



(40) 



In the following we restrict 7 to be a positively oriented circumference going through ±i and parameterize it as 
z(9) = —a + vl + a 2 e t6 , with 6* S [0, 27r) and a > a constant. Changing variables according to the Mobius 
transformation s(t) = and denoting Qi(s) = A + Bs 2 , Q2(s) = C + Ds 2 we can write (|4(JI) in the form 



7Ti?l 



sgn!ft 



s-l 



i(s+l) 



( S -l)4 



1 J S ( 7 ) 



( S -l) 



(s-l) 4 



^1 A( T ) VWJOaW 



3(1) 



(41) 



where 5(7) -the image of the circumference 7- is a straight line through the origin with slope — ~ < 0. The sign in 
the numerator in the last integral, for each s G C, is shown in fig. [IJin the different regions in the (p, r) parameter 
space. We now proceed to compute the integral I appearing in (|41|l for the regions IA, IB, IIA, IIB, and III, and for 
the hyperbola r 2 = p 2 + 1. 
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Region IA: Here A < 0, B > 0, C > 0, D < 0, and the sign in the integrand is positive. We choose a — > and 
parameterize s(7) as s = cr(i — e) [er 6 M, e — * + ]. We have 



/=-• 



da 



nRi J R y/(A- Ba 2 ){C - Da 2 ) - ie 

with e — > + . Recalling our choice of branch for the square root, changing variables according to u = *J — -jc, and 
noticing that < < 1, this integral can be written as 



1 1 f du 2 1 / / 4p . 



«RiVBCh / {l + u 2 ){l+ AD u2) nR lVi l +p) 2^ T 2 ^(l + p) 2 -r 



Region IB: Here A > 0, B > 0, C > 0, D > 0, and the sign in the integrand is positive. We choose a — ► oo and 
parameterize s(7) as s = cr(— 1 + ie) [er e R, e — > + ]. Then 



nRi Jr ^/(A + Ba 2 )(C + Da 2 )-ie 



with e — > + . After the change of variables u = w ^<t (and since < < 1) this integral becomes exactly <|42[) . 

Region IIA: Here A < 0, -B > 0, C > 0, D > 0, and the sign in the integrand is negative. Letting a — ► and 
parameterizing s(7) as s = er(i — e) [a e M, e — > + ] we obtain (with e — > + ) 

rftr 



This last integral can be split in two pieces: 



We find then (see formulas 3.152-3 and 3.152-6 of |13j 1 



1 jff./q+y-^U'jfC K-'r"'i. (43) 



Region IIB: Here A < 0, i? > 0, C > 0, D > 0, and the sign in the integrand is positive. With the choice a — > oo 
and the parametrization s = <r(— 1 + ie) [cr G R] for s(7) we arrive at 

1 /" do- 



TT-Rl 7r v /(A + -B(7 2 )(C + L'cr 2 ) - ie 

with e — ► + . This can be written as 

_ 2i /•/"¥ dcr 2 f 00 r/o- 



tt-Ri Jo ^/-{A + Ba 2 ){C + Da 2 ) ttRi J ^/(A + Ba 2 )(C + Da 2 ) ' 
which gives again (|43() . 

Boundary between regions IIA and IIB: This is the hyperbola r 2 = 1 + p 2 . Parameterizing s(^) as s — a(i — e) 
[a e E, e -> 0+] we get 

i f da 



ttRi Jw y/(A- Ba 2 )(C - Da 2 ) + ie 
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(with e — > + ) which reduces to the previous two cases. Substituting r 2 = 1 + p 2 we thus obtain 

2 



TrR ly fp \V2J A-r:R ly fWp 



1 + i 



Region III: Here A < 0, £> > 0, C < 0, D > 0, and the sign in the integrand is negative. Letting a 
parameterizing s(7) as s = a(i — e) [a G E, e — > + ] we obtain (with e — > + ) 



da 



With the change of variables u = \l — %a and noticing that < -jjj < 1, this becomes 



and 



I = 



du 



2i 



*Ri V^AD U / {1+u 2 ){1+ bc u2) kRi ^Jt 2 - (1 - P) 



-.K 



Ap 



VIII. APPENDIX IV 

We want to compute the double integral 137(1 . In order to do this we first integrate in t\. Fixing t% we find that 
the poles coincide with those of 1(36(1 . The arguments showing whether they are inside or outside 7 are the same as in 
Appendix III. The residues are now 

t+ , t 2 [p 2 (t 2 -l) 2 +A l prt 2 (t 2 ~l)^2(l + 2r 2 )t 2 } 



1 [4t| + (-p + pt\ + 2iTt 2 ) 2 WU'i + (-p + pt'i + 2irt 2 ) 2 ' 
r : t 2 [p 2 (t 2 -l) 2 +4zpTt 2 (t 2 -l)-2(l + 2T 2 )t 2 ] 

1 [At 2 + (-p + pt\ + 2iTt 2 ) 2 ] y/At 2 + (-p + pt 2 + 2iTt 2 ) 2 ' 

Integrating in t\, ((37(1 becomes then 



3/ 4T I d t S -^M^L + El + P f2 + 2zT ^ 2 J [^^ 2 ~ ^ + = 1} ~ 2(1 + 2T2)t2] ^ 'f (4 1. ) 

°\ 7ri?iA/ 7 [4t 2 + (-p + pt 2 + 2irt) 2 ] 2 yjAt 2 + (-p + pt 2 + 2irt) 2 

The same Mobius transformation used above allows us to write ((44(1 as 



^ 1 — n — rTT / ds sgn 3? 



27Ti?lA 



»(7) 



i(s 2 - 1) 



(.s 2 - l) 2 [2p 2 (l + s 2 ) 2 + 4pr(s 4 -!) + (! + 2r 2 )(s 2 - l) 2 ] 
Ql^)Ql{s)^/Qi{s)Q 2 {s) 



We compute it only in regions IA and IB, remembering that the sign in the integrand is that shown in fig. [5J 

Region IA: Here A < 0, B > 0, C > 0, D < 0, and the sign is positive. We choose a — > and parameterize 5(7) as 
s = a(i — e) [a G K, e — > + ] getting 

J (1 4- a 2 ) 2 [2p 2 (l - a 2 ) 2 + 4pr(a 4 - 1) 4- (1 + 2r 2 )(a 2 + l) 2 ] 
da- 



2irR 1 \J R (A ~ Ba 2 ) 2 {C - Da 2 ) 2 ^-{A - Ba 2 )(C - Da 2 ) 

The integrand consists of a rational function of a divided by the square root of a fourth degree polynomial and the 
integration extends over the real axis; hence, it can be written in terms of complete elliptic integrals. The way to 
proceed now is to perform a partial fraction decomposition of the rational part and write it as a sum of integrals of 
the four following types 

dx f dx f dx f dx 



sj(a 2 + x 2 )(b 2 + x 2 ) 3 Jr y/(a 2 + x 2 ) 3 (b 2 + x 2 ) Jr y/(a 2 + x 2 )(b 2 + x 2 ) 5 Jr y/(a 2 +x 2 ) 5 (b 2 + x 2 ) 

These integrals can be found in ^| (see equations 3.158-2, 3.158-4, 3.162-4, and 3.162-2 in that reference). The result 
is the 1/A contribution appearing in equation ((13(1 . 
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Region IB: Here A > 0, B > 0, C > 0, D > 0, and the sign in the integrand is positive. With the choice a — > and 
the parametrization s = a(—l + ie) [a G M, e — > + ] we obtain 

cJ g /• ^ - l) 2 [2p 2 (l + ^ 2 ) 2 + 4pr(a 4 -!) + (! + 2r 2 )(q 2 - l) 2 ] 
^1 2ttJ2iA 7 R CT (A + Ba 2 ) 2 {C + Da 2 ) 2 ^{A + Ba 2 )(C + Da 2 ) 

It can be shown that this integral gives again the result found in region IA. 
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